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1^ ' Abstract. We study the quantum analogs of tops on Lie algebras so(4) 

' and e(3) represented by differential operators. 

o 

m 

^ ' 1 Introduction 

I The commutative rings of differential operators play an important role in 

mathematical physics. The problem of their description, in the case of one 
independent variable, was set up and solved in the works by Schur [1] and 
Burchnall-Chaundy [2]. Later on, these results have found numerous ap- 
J> ' plications in the theory of integrable equations of the Korteweg~de Vries 

^ ■ type; we mention only the description of finite-gap operators obtained by 

I Novikov, Krichever and others [3, 4], and the algorithm of checking necessary 

' integrability conditions developed by Shabat and others [5]. 

Q\ • In the case of two independent variable, the generalizations in different 

, directions are possible. In general, the problem on a pair of commuting 

operators [H, K] = is not very well posed and the question about addi- 
tional assumptions arises. Apart from the conditions suggested by the logic 
of the problem itself (for instance, commutativity of a triple of operators), 
^ I it is also natural to select the classes of operators which are important in 

' applications. Here we consider one such class related to the quantum tops 

on so(4) and e(3). The corresponding operator if is of a very special form 

H = a{x)Dl + 2/(x, y)D.,Dy + b{y)Dl + . . . 

where a, h are fourth degree polynomials, /(x, y) is a biquadratic polynomial 
and dots denote the lower order terms, also with special coefficients. This 
form appears if we consider a Hamiltonian which is quadratic with respect 
to the generators of Lie algebra, for instance 

H = ({/, AfJ) + (T/, BV) + 2(C/, FV) + (c, {/) + {d, V) 
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in the case of so(4), and use a differential representation such that generators 
U, V are replaced with differential operators of first order with quadratic 
coefficients (see equations (6), (8) below). The order of operator K does 
not exceed 4 in the known examples and, respectively, the degree of its 
coefficients does not exceed 8. The original classical version of the problem 
deals with Hamiltonians polynomial in momenta which are in involution 
with respect to Poisson-Darboux bracket, instead of commuting differential 
operators. The conditions of "commutativity of the principal part" [6], that 
is, cancellation of the leading terms in the commutator, coincide for the 
classical and quantum cases. In more details, this relation is illustrated by 
the diagram 

f , quantization on the Lie algebra r -y y-, 

\H, K I = (J > [H , K\ = L) 



R* 



R* 



r 1 n quantization in the Darboux coordinates r " ^ n 

{Hn, Kd} = > [Hd, Kd] = 

Here R denotes a representation of the Lie-Poisson bracket in the Darboux 
coordinates, and Hd is the image of the Hamiltonian H under the pullback 
map, that is Hd = R*{H) = HoR. Analogously, R denotes a representation 
of the Lie algebra of differential operators and Hd = R*{H) = H o R. The 
horizontal arrows correspond to the passage from the Lie-Poisson bracket to 
the commutator on the Lie algebra and from the Poisson-Darboux bracket 
to the Heisenberg algebra 

This passage, on both levels, is performed by choosing a suitable ordering 
of the monomials of Hamiltonians H, K, which is equivalent to adding the 
lower order terms ("quantum corrections"). 

Let us remind that the description of the classical tops admitting an 
additional first integral is still an open problem. In the ,80(4) case, the 
necessary condition were derived by Veselov [7] (see [8] for a more precise 
result) which define a certain varieties in the parameter space. However, 
the number of parameters in the known examples is less at least by 1 than 
the dimension of these varieties. Some classification results are presented in 
paper [9]. Is it possible to make progress in this problem, or, more precisely, 
in its quantum version, by use of the language of differential operators? As 
a first step, in this paper we construct the differential operators which cor- 
respond to several most important examples. In particular, these examples 
suggest that probably the simplest necessary condition of integrability is 
that the discriminant of the leading part of H is factorable: 

f{x, yf - a{x)h{y) = w{x, y)w{x, y) (1) 
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where w, w are, as well as /, biquadratic polynomials. Notice that this 
relation is important in the problem of separation of variables [10]. It is not 
clear now, how to prove this property in the general setting, using only the 
assumption that there exists a commuting operator K of arbitrary order. 
From the computational point of view, the advantages of the differential 
representation is that it takes care of such issues as ordering in the operator 
algebra and Casimir operators. Notice that representation R is defined on 
a fixed level set of Casimir operators and their values occur in H£,,Ko as 
independent parameters (their number is equal to the rank of Lie algebra, 
that is 2 in so(4) and e(3) cases). Splitting of the commutation relation with 
respect to these parameters makes its analysis more easy. 

Contents of the paper are the following. The so(4) case considered in sec- 
tion 2 includes the tops of: Schottky-Manakov [11, 12], Stekloff [13], Adler- 
van Moerbeke [14, 8] and Sokolov [15]. The quantization of the Schottky- 
Manakov top was studied in paper [16] . The cases of Adler-van Moerbeke 
and Sokolov are more complicated (the operator K is of fourth order) and, 
up to our knowledge, their quantization is obtained here for the first time. 

In section 3, we consider tops on e(3): the Clebsch top [17], the Kowalevskaya 
top [18] (quantization is due to Laporte [19]), the Kowalevskaya gyrostat and 
the Goryachev-Chaplygin case. The quantization of models on e(3) is rather 
well studied by Komarov and others [20, 21, 22, 23, 24]. 

Certainly, the above examples do not exhaust all known cases. An ac- 
count of the known results and a detailed bibliography can be found in book 
[25]. Other directions in the problem on commuting differential operators 
were developed by many authors, we mention only papers [26, 27]. 

Notations. We denote x = x\, y = X2 for short. The definitions of Poisson 
brackets and commutators do no include vanishing brackets. The subscripts 
z, J, k always denote an arbitrary permutation of 1, 2, 3. The imaginary unit 
is denoted as i. The notations 

a = (01,02,03), {o.,h) = 0161 -I- 0262 + 03^3 

are used. The following polynomials are often used in writing a compact 
form of Hamiltonians in Darboux coordinates: 

W{a, b, c; x, y) = o(x^ - l)(y^ - 1) - b{x'^ + l){y^ + 1) + 4cxy, 
R{a, b, c; x) = W{a, b, c, x, x) = {a - b){x^ + 1) + 2(2c -a- b)x'^ . (2) 

Notice that operator VF(o, 6, c; x, y)DxDy is form invariant under the inver- 
sion X— )-l/x, y^l/y. 
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2 Tops on so(4) ~ so(3) so(3) 

2.1 Darboux coordinates and operator representation 

The Lie-Poisson bracket on so{3) and its Casimir function are defined by 
equations 

{U„Uj} = ie,jkUk, {U,U) = sl (3) 

where U denotes the vector {Ui, 1/2,113). A representation in Darboux co- 
ordinates is given by relations 

Ui = -^{x'^ - + Six, U2 = -^{^'^ + +isix, 

U3 = -xpi + Si. (4) 

In order to quantize, we pass to the operator algebra 

[Ui, U,] = i EijkUk, {U, U) = jiUi + 1) (5) 

which is represented by diff'erential operators as follows: 

Ui = -^(x^ - l)D^+jix, U2 = -^(x^ + 1)D^ + ijix, 

U3 = -xD^+ji. (6) 

The Lie algebra so(4) ~ so(3) © so(3) is represented as the direct sum with 
a double under the following renaming: 

U^V, x^y, pi^P2, si^S2, U^V, Dx^Dy, ji^j2- 

The tops on so(4) are defined by Hamiltonians of the general form 

H = (U, AU) + {V, BV) + (U, FV) + (c, U) + {d, V). (7) 

The matrices A and B can be chosen diagonal without loss of generality. 
The use of representation (6) yields a differential operator of the special 
form 

Hd = a{x)Dl + f{x, y)D,Dy + h{y)Dl 

^"'^ ^ a'{x) + 2c(x) + i2fy{x, y)) 



^^^6'(y) + 2d{y)+hUx,y))Dy 

o b 

+ 2jic'{x) + 2j2d'{y) + hhfxyix, y) + k (8) 
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where k is an arbitrary constant and 

a{x) = aoix^ + 1) + a2X^, h{y) = bo{y^ + 1) + 62y^ 
c(x) = C2X^ + cix + Co, d{y) = d2y'^ + diy + do, 
f{x,y) = f22X^y'^ H h /do- 

The correspondence between (8) and (7) is given by relations 

A = diag(2ao, -2ao, 02) + al, B = diag(26o, -26o, ^'2) + 1^1, 
c = 2(c2 - Co, - i(c2 + Co), ci), d = 2(^2 - do, - 1(^2 + 4), di), 

(/OO — /02 - /20 + /22 i(/00 + /02 — /20 — /22 /2I — /oi 

i(/oo - /02 + /20 — 722) — /oo — /02 - /20 - /22 — i(/2i + hi) 
fl2 — fw — i(/l2 + /lo) /ll 

3k = (3a + a2)ji(ji + 1) + (3/3 + 62)72(72 + 1) 

with arbitrary constants a, /3. 

The fohowing remark is useful for computing spectra. It is about the re- 
lation between differential representation (6) and the matrix representation 
[28, p. 88] 

Ui±iU2 = U±, U3\mJ)=m\m,j), 



U+\m,j) = A/(i-m)(i + m+l)|m+l,j), 



U-\m,j) = A/(i + m) (j - m + 1) |m - (9) 

In the physical case we consider an orbit j = const such that j and m are 
simultaneously integer or half-integer and m takes the values from —j to j, 
which corresponds to the matrices of the finite size {2j -|- 1) x {2j + 1). For- 
mally, one can waive these restrictions and consider the infinite tridiagonal 
matrices which are interpreted as the difference operators 

a{m)Tm + b{m) + c{m)T~^ 

where : m 1— t- m -|- 1 is the shift operator (the matrix with units on the 
subdiagonal). A conjugation A — f~^Af allows to get rid of square roots, 
and the representation of so(3) by shift operators appears: 

Ui = -\{rn - j){m + j + 1)T„ + 

C/2 = - i(^(m - j){m + j + l)Tm. + r~i) , 

U3 = m. (10) 

Another difference representation can be obtained directly from (6) via the 
change ("Fourier transformation") 

mTm, X 
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which preserves the Heisenberg algebra: 

[D^,x] = l, [mTm,T~^] = l, 
and yields the operators 

( m\ „ 1 

C/i = yT„,+ (j + l--jr^i, 

c/2 = -iyr^ + i(i + i-y)T-S 

C>3 = j - m + 1. (11) 

One can check that representations (10) and (11) are equivalent up to a 
certain linear transformation from 5*0(3) and the conjugation by operator 
g{7n)T~^~^ where function g satisfies relation g[m + 2)/g{m + 1) = 2j — m. 

2.2 Schottky— Manakov top 

The commuting operators are of the form 

H = -alul - alul - apl - ajV,^ - a^V^ - alvf 

+ 2a2a3UiVi + 2a3aiU2V2 + 20102^^3^3, (12) 
K = aiUiVi + a2U2V2 + a^tj^%. (13) 

These can be obtained from the classical Hamiltonians in involution by the 
simple exchange U ^ U , V ^ V ^ because the problem of normal order- 
ing does not appear in this example (each monomial contains commuting 
variables only). In contrast to some other examples (Stekloff, Kowalevskaya 
tops), no linear terms can be added to the Hamiltonians, both in classical 
and quantum cases, as one can prove by a direct computation with indeter- 
minate coefficients. Applying map (6) yields the following pair of commuting 
differential operators: 

= r{x)Dl + 2zD^Dy + r{y)Dl 

- [^^r'{x) + 2j2Zy^ D.^ - (^^r'(y) + 2jiz..) Dy 

- ^(ii(ji + 1) +i2(j2 + l))(a? +al + ai), (14) 

AKd = wD^Dy - j2WyDx - jlWxDy + ilj2Wxy (15) 

where, using the notations (2), 

r{x) = -R{a\,al,al;x), 

z = z{x, y) = ^^(0203, asai, 0102; x, y), 

w = w{x, y) = W{ai,a2, 03; x, y). 
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The constant in the last line of (14) (k in equation (8)) has no effect on the 
commutativity and determines the energy shift of the ground state of the 
system. 

Alternatively, operators (14), (15) can be obtained by another route 
along the diagram from Introduction. First, one should pass to the Darboux 
coordinates in the classical Hamiltonians H, K, accordingly to formulae (4). 
This yields 

AHd = r{x)p\ + 2zpip2 + r{y)pl 

- (sir'(x) + 2s2Zy)pi - {s2r'{y) + 2siz^)p2 

+ 2siS2Zxy + |r"(x) + |r"(y) -\{sl + sl){al + al + oi), (16) 

AKd = WpiP2 - S2WyPi - SiWxP2 + SiS2Wxy (17) 

where r, z, w are the same as above. Next, pi are replaced by in quadratic 
terms and the lower order terms are taken with indeterminate coefficients 
(preserving the degree with respect to x^y). So, one searches for the com- 
muting differential operators of the form 

Hd = r{x)Dl + 2z{x, y)DxDy + r{y)D^ + a{x, y)D^ + h{x, y)Dy + c{x, y), 
Kd = w{x, y)DxDy + /(I, y)Dx + g{x, y)Dy + /i(x, y) 

where r, z, w are given, a, 6, c, /, h are polynomials with indeterminate co- 
efficients, and numbers above arguments show the corresponding degrees. 
The resulting system of equations is rather bulky, but its solution (14), (15) 
can be easily recovered by any system of computer algebra. Thus, both 
methods of quantization turn out to be equivalent, as one should expect. 
In the situation, when H is given and K is to be found, each method has 
its own strong and weak points from the point of view of computation com- 
plexity. The use of the Lie algebra generators gives a simpler equations 
for the coefficients, but its implementation is more difficult because it deals 
with noncommutative variables. The language of differential representation 
is more flexible and universal, but the resulting system for the coefficients 
of the operators is slightly more cumbersome. 

2.3 SteklofF top 

The classical Hamiltonians are of the form (a = aia2Cfi) 

H = -o? {\ul + + ^C/| ) + 2a(ai[/iFi + 02^73^2 + asf^s^s), 

\a{ al ) 

K = 2a (—UiVi + —U2V2 + —U;V^ - alV^ - alV^ - ajvi. 
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Like in the previous example, the problem of ordering does not appear and 
the commuting quantum Hamiltonians are obtained just by placing hats 
over U, V. It is easy to check that the following linear terms can be added 
to the Hamiltonians, both in the classical and quantum cases: 

H H - (iia2a^Ui - /32aia3f^2 - /33aia2C^3, 
K ^K + PiVi+hV2 + PzVi, 

with arbitrary /3j. However, we assume /3j = for the sake of simplicity. 
The passage to the differential operators yields 



AHd = nDl + 2zD^Dy - \^^^i±^r[ + 2j2Zyj - 2jiz^Dy 

I o ■ ■ I ~ 1) // ^ ■ r ■ I i\/ 2 2 , 2 2 I 2 2\ 

+ ^Jij2Zxy H ^ - gJiUi + l)(aia2 + a2"3 + "3«i) 

AKd = 2wDxDy + r2D^y - 2j2WyDx - (^^^^^^2 + 2ii^^x) Dy 

J2(2i2 - , 2 , 2 , 2n 



2ii - 1 , 



+ 2jij2Wxy H r2 - -J2(j2 + + ^2 + "s) 



where 



/\ 0/22 22 22 \ 

n = ri{x) = -R{a2a^,a^ai,aia2;x), 
r2 = r2{y) = -R{al,al,al; y), 
z = z{x,y) = aia2a3W{ai,a2,a3;x,y), 
w = w{x, y) = W{a2a3, asai, aia2; x, y). 

The same operators are obtained via the quantization in the Darboux vari- 
ables, starting from the Hamiltonians 

AHd = ripj + 2zpip2 - isir[ + 2s2Zy)pi - 2siZxP2 

g2 4 

+ 2siS2Zxy + -jr'l - -sjiajal + a^al + alaj), 
AKd = 2wpip2 + r2p\ - 2s2WyPi - {s2r2 + 2siWx)p2 
+ 2siS2Wxy + -jr2 - -sliaj + al + aj). 
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2.4 M. Adler— van Moerbeke top 

The classical Hamiltonians are of the following form (the parameters are 
related by the constraint Ai + A2 + A3 = 0): 



H = ^(-9A|A^.C/:f + 6A,Afc(A, - Xi){\k - Xi)UiVi 
1=1 

+ A,Afc(4A2-A,Afc)F,2), (18) 
K = ^Y. - >^3)U^y^Vj + J^(Ai - A,)(Ai - Xk)U^Vl 

- 9 E E ^^^'^^^^'^ + ^ E E(3A^ - - ^i)^^ (19) 

I i I i 

The quantum Hamiltonian H is obtained from H by fixing hats. However, 
K contains monomials with noncommutative generators and the problem 
of ordering arises. The result of direct computation with indeterminate 
coefficients is the following expression for K: 

^ = E - Aj)(c>iV-y/ + UiVjV^Vj + u^vfVi) 

+ E(^^ - - Afc)^>^V^^' - 9(E + ^) E AjAfef/iV- 

i I i 

+ ^E^'E(3\^-\'-^i)^^ (20) 

I i 

Thus, the quantum Hamiltonian differs from the classical one by ordering in 
the first sum and a quantum correction in the third sum (recall, that 
is the Casimir function). 

The following expression for Hd is obtained by passage to differential 
operators according to general formula (8): 

iHD = ri{x)Dl + 2zD,Dy + r2{y)Dl 



^(9ii(ji + 1) + J2(j2 + 1))(A? + A1A2 + xlf (21) 



3 

where 

ri{x) = -9R{XlXlXlXlXlXlx), 

r2{x) = \ri{x) + 4AiA2A3i2(Ai, A2, A3, x), 

z = z{x,y) = 3W{ni,i22,i^3,x,y), m = {Xj - A^)(Aj^ - A|). 
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The structure of operator is rather simple: 

AKd = g^D^Dl + ciggyD^Dl + ca^^^-D^ 

+ {c?.99yy + CA9l)DxDy + {c;,gygx + ceggxy)Dy 

+ {ciggyyy + C8gygyy)Dx + {cggyyg^ + c^gygxy + cuggxyy)Dy 

~l~ (^129yyy9x ~\~ C.1^9yy9xy ~\~ ClA9y9xyy ~\~ Cl599xyyy ^ 

but the coefficients are cumbersome. Here g denotes the polynomial 

5 = (Ai - \2){xy^ + 1) + 3A3(x + y)y 

related to the coefficients of Hd by relation 

- ri{x)r2{y) = 36AiA2A3(Ai - A2)(A2 - A3)(A3 - Ai)^^, 

where g = (Ai — A2)x(2; + y) — Xs{x'^y + !)• The coefficients Cj are, in turn, 
polynomials in the parameters ji,j2- 

ci = -2(i2 - 1), C2 = -2ji, 

C3 = ^(3jl(jl + 1) + J2(j2 - 3) + 2), C4 = -jlUl + 1) + i2(j2 " 1), 
C5 = 2ii ijl + J2), C6 = -2ji {ji - J2 + 2), 

C7 = (9jl(jl + 1) - J2U2 + 3) + 4), C8 = + 1) - J2(J2 - 1)), 

eg = - y (3ji + 4ji j2 + j'i + ji + j2 ) , Clo = 2ji {ji {ji + 1 ) - J2 ( j2 - 1 ) ) , 

cii = -f - 4jii2 + jl + 5ji - 7j2 + 4), 

ci2 = ^(9j? + 3jii2 - j| + 6ii + 1), 

ci3 = ^(-3j? + 5iij2 + 3j| - 4ji + 2i2 - 1), 

ci4 = ^(-3j? - 5iij2 + 3i| - 2ji - 8i2 + 1), 
b 

ci5 = ^(9jf - 3jii2 - j| + 12ji - 6j2 + 7). 
2.5 Sokolov top 

In this example, it is convenient to use the variables 

mi = Ui + Vi, ni = Ui- Vi. 
The Sokolov top on so(4) is defined by the Hamiltonian 

H = \m\ + \ml + m| + m3(ani + ^n2) -\{o? ^ f)n\ (22) 
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and additional fourth order integral is of the form 

K = ml[2H - m| + (/3mi - am2)^ + [ani + /3n2)^) 

+2m-i[ami + /3m2 - (a^ + /3^)n3) (mini + m2n2). (23) 

The quantum Hamiltonian is obtained by symmetrization of noncommuta- 
tive monomials: 

H = hhl + + ml + [m3,ani + /?n2]+ - ^(a^ + /3^)n^ (24) 

where 

[a, 6]+ = -(a6 + 6a). 

The operator K is obtained by symmetrization as well, with some special 
weight coefficients. It should be noted that there exist homogeneous poly- 
nomials in the algebra so(4) which vanish identically in virtue of the com- 
mutation relations, so that the form of the operator K is not unique and 
one can try to simplify it by adding such polynomials. We write down one 
of possible versions. Let us denote 

^a,fe,c("i, /) = ampm + h[w? , /^]+ + c[/, mfm]^ 
+ {l-a-b-c)fn?f, 
B{m, /, g) = mfg + gfm, 

then the operator 

K = ^1 ,o,| ("^3, iTT-i) + ^i,o,f ("^3, "^2) - («^ + 0^)'rh\ni 

+ Ai 1 , (ms, /3mi — 07712) + A_a ^ 7 (ms, ms + ahi + /3n2) 

2 ' 4 ' 4 ' ' 2 

+ B{m^, ami + /37n2 - (a^ + /3'^)hs, rhihi + m2"'2), 

- ^["13%, ["i3,ami + /3m2]] (25) 

commutes with and coincides with K if all variables are commutative. 

In this example, it is also convenient to modify slightly our differential 
representation of so(4), by changing signs y —y, Dy — ?• —Dy. This brings 
to the operators which are symmetric with respect to y. Let us denote 
(in general, the parameters a, (3 are complex) 

^i = a + i/3, ^2 = -a + i/3, r{x) = x{iix + l){x + ^2), 
z = z{x, y) = iixy{x + y) + {x + yf + 2(1 - 66)a;y + 6(a; + y) 
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then H is written in the form (8) as follows: 
2Hd = r{x)Dl + zD^Dy + r(y)L>2 



+ luiUi + 1) + J2(j2 + i))(eie2 + 4). 

Notice that the polynomial r is of degree 3, because Hamiltonian (22) takes 
the form (7) with nondiagonal matrices A, B if one returns to the variables 
U, V . The second operator is too bulky and we present explicitly only the 
leading terms: 

Kd = w^{xD^ + yDyfD^Dy - 2j2X^wWyDl - 2jiy'^ww^D^ 

- 2xw{{ji - l)xwx + (ji +j2 - l)w + (2j2 - l)yWy)DlDy 

- '2yw{{2ji - l)xwx + (ji + J2 - l)w + (i2 - l)ywy)DlDy + ... 

where w = £,ixy + x + y + ^2- 
2.6 The classical limit 

The Planck constant is introduced by simple scaling of the generators, so 
that commutation relations (5) are replaced with 

\Ui, Uj] = i hei.kUk, {u, u) = h^h {ji + 1) 

and representation (6) is replaced with 

Ui = - ^)Dx + jix) , U2 = ^{-^{^^ + + i Ji^;) , 

Us = h{-xD^+ji). 

The equations for the variables V are changed analogously. The passage to 
the classical limit for any quantum operator A is defined according to the 
formula 

A = lim e"i(Pi^+P^^) ( A e^(Pi^+P22^)') I . s, , 

in particular the commutator bracket and the Casimir functions for each 
copy of so(3) are mapped into the Lie-Poisson bracket (3), and the formulae 
for the generators themselves are mapped into the Darboux coordinates 
representation (4). 

Applying of this procedure to Hd, Ko gives the same expressions for Ho, 
Ku as the intermediate passage to the Darboux coordinates in the classical 
Hamiltonians H, K. This is guaranteed by the "correspondence principle" 
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which is invariant, that is it does not depend expUcitly on the choice of rep- 
resentation of the algebra so(4). The check of the correspondence principle 
is trivial for the Schottky-Manakov and Stekloff tops, because there are no 
quantum corrections in these systems. In the Adler-van Moerbeke case, the 
correct passage to the classical limit is achieved by changing one term in the 
expression for K (20): 

I i I i 

After this, all terms become homogeneous with respect to h, and the limit 
/i — 7- gives rise to the classical Hamiltonian (19). 

It should be noted that the Casimir functions + 1) are of the 

quantum nature, because ji take integer/half-integer values. The passage to 
the limit ^ brings to the classical (finite) quantities Sj = hji. On the 
other hand, if we consider spins then the values ji are finite and therefore 
Si — > 0, in accordance with a statement that the spin is a pure quantum 
concept. 

3 Tops on e(3) 

3.1 Darboux coordinates and operator representations 

The Lie-Poisson bracket on e(3) is of the form 

{Mi,Mj] = -EijkMk, {Mi,-fj} = -Eijklk, {lulj} = 

and the Casimir functions are 

(M,7) = Z, (7,7) =0^ 

where M = (Mi,M2,M3), 7 = (71,72,73)- We use the following represen- 
tation in the Darboux coordinates: 



Ml 


= - - 


^(y' - 1)P2 




M2 


= -\{x' + l)pi - 


\{y' + 1)P2 


-i^(x-y) 


M3 


= \{xpi + yp2), 






71 


1 — xy 
= a , 72 = 
x-y 


. 1 + a^y 
1 a , 

x-y 


x + y 

73 = a 

x-y 



(26) 
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A representation with real Darboux coordinates [29] should be mentioned 
as weh: 

Ml = -piqiq2 + -P2{qi - q2 - I) + ^ . 2 ^ 2^ ' 
2 2a{qf + q^) 

M2 = P2qiq2 + T,Pmi - 92 + 1 + ^ . 2 , 2\ ' 

M3 =pii?2 -P2Q'i, 
^ 2agi ^ 2ag2 ^ +^2-1) ^r,7x 

('Zf + 'zi + l)' {<ll + 4 + ^y {ll+<ll + l)' ^ ' 

Quantization replaces the Lie-Poisson bracket with the e(3) commutator 
[Mi, Mj] = i EijkMk, [Mi,7j] = i Sijklk, = (28) 

and the Casimir operators are 

(7,M) = /, (7,7) = a'- (29) 
This operator algebra admits the following representation: 

Ml = \{l - x')D, + 1(1 - y')Dy + ^{x - y), 

M2 = ^(1 + x')D, + ^(1 + y')Dy -il-{x-y), 

Ms = xD^ + yDy, 

^-xy ^ . 1 + xy „ x + y 

71= a , 72 = ia , 73 = . (30) 

X — y X — y X — y 

Notice that all operators are invariant with respect to the change 

X y, o — > —a. 

A matrix representation can be obtained by introducing the basis func- 
tion 1^) = \m,n) = (x + y)"^{x — 7/)"". This choice is motivated by the 
denominator of generators 7^ in representation (30) and the symmetry ar- 
guments. Easy computation yields {M± = Mi it iM2, 7± = 71 ±172) 

Msim, n) = {m + ?i)|m, n), M_|m, n) = 2m\m + 1, n — 1), 
M_|_|m, n) = — (n + — j |m + l,n) — —\m— l,n + 2) H — |m,n + 1), 
73|m,n) = a\m + l,n — 1), 7„|m,n) = 2a\m,n — 1) , 
7+|m,n) = — \m,n+ 1) — — \m + 2,n — 1). 
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In order to pass to the classical limit, the Planck constant is introduced 
as follows: 

[Mi, Mj] = i hsijkMk, [Mi,7j] = i heijk^fk, = 0, 

and the operator representation is replaced by equations 

Ml = - x^)D., + ^(1 - y^)Dy') +^ix- y), 

M2 = + x^)D., + ^(1 + y')Dy') -il-(x-y), 

M3 = h{xD^ + yDy), 

1 — xy , 1 + xy X + y 
11 = a , 72 = ia , 73 = a . 

X — y X — y X — y 

The formula for the classical limit is analogous to the so(4) case: 
A = lim e"i(Pi^+P2?/) ( ^ f^i{pix+P2y)\ 

however, notice that here the Casimir operators (29) are pure classical. Ap- 
plying this procedure to the generators Mj yields the bracket (26). In the 
Kowalevskaya case this procedure results in changing of a coefficient in op- 
erator k (35) (cf [19]): 

k = ^{k+k- + kA+) + Ah^{Ml + M|). 
3.2 The Clebsch top 

There is no problem of ordering in this case and the quantum top is defined 
by the Hamiltonians 

H = -Y,{^f + \af), (31) 

i=l 

^=^E(^*^^^-^|J)' ^ = ^i^2A3. (32) 
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The use of representation (30) yields the following commuting differential 
operators: 

2Hd = -{x - yfD^Dy + ^(x - y){D, + Dy) 

+ . " + a^(Ai + A2 + A3), 
[x-yY 

SKd = r{x)Dl + 2zD^,Dy + r{y)Dl 

+ (^/(.) + Lzy)D. + {^r'iy) - Lz.)Dy 

+ (Ai - A2)^(x - ,)2 - (A, - A,) - y') + .^ K^^My)^-^^ 

where 

r(x) = i?(Ai, A2, A3; x), z = z{x, y) = W'(Ai, A2, A3; x, y). 
Notice, that in this case identity (1) takes the form 

z" - r{x)r{y) = 4(x - yf {\i\2 + A2A3 + A3A1) 

+ 4(x - y)2p^(A2A3, A3A1, A1A2; X, y), 

therefore the last term in Kd partially cancels. The form of classical Hamil- 
tonians in the Darboux coordinates is analogous, with slightly different co- 
efficients: 

21 

2Hd = -(x - yfpiP2 H {x- y){pi +P2) 

a 

+ " ^,0 + a^(Ai + A2 + A3), 
(x - yy 

8Kd = r{x)pl + 2zpiP2 + r{y)pl 

- -(r'(x) - 2zy)pi + -{r'{y) - 2z^)p2 
a a 

+ 4(Ai-A.)4(.-,)' + a^ '-'ffa'-/' . 

(x — y)^ 

3.3 Kowalevskaya top 

The classical top is defined by Hamiltonians 

H = ^ {Mf + M| + 2M|) - , K = k+k^ (33) 

where 

k± = (Ml ± iM2)2 + 71 ± 172- 
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The expressions in the Darboux coordinates (26) are: 

2H = x^pI + {Axy - x^ - y'^)piP2 + y'^pl + — (x - y){pi +P2) + a— — 

a X — y 

( f 2 2 2Z \^ 2axy\ ( , ,9 2a 

K=\\xpx^yp2 \x-y)\ (Pi+?'2) ' 



a J X — y I \ X — y 

The quantum version of the Kowalevskaya top is of the form [19, 22] 

H = \{Ml + M| + 2M|) - ^71, (34) 
k =]^{h+k- + k-k+) + 4.{Ml+M^) (35) 

where 

k± = (Ml ±iM2)^ + 7i ±172- 

The use of representation (30) yields the following commuting differential 
operators: 

2Hd = x^Dl + (4xy - x^ - y^)D^Dy + y^Dl, 

+ -((a + l)x - ly)D^ + -{lx + {a- l)y)Dy + a— , 

a a X — y 

kD=[f-^,g' + ^V+ 4[/, gfg]^ - 2[fg, gf]^ -2[f, gY 
'- X — y X — y' 

where [a, ^J"*" = ^{ab + ha) and 

f = x'^D^ + y^Dy-J-{x-y), g = D^ + Dy. 

It is worth noticing that Hamiltonians (34), (35) admit the following gener- 
alization (Kowalevskaya gyrostat) [23]: 

2H = Ml + M| + 2M| - 71 + cMs, 

k = ]^{k+k^ + kSk+) + ^{mI + M|) - 2c{mI + M|)M3 

+ 2c^M| + c(c^ + 1)M3 - 2cMi73 - c^7i - i C72. 

3.4 Goryachev— Chap ly gin case 

The quantization was considered in [21]. The Hamiltonians 

H = Ml + M| + 4M| - 71 + cMs, 

k = 4(Mf + M|)M3 + 2Mi73 - 4cM| + (1 - c^)lvh + C71 +172 
satisfy the relation 

[H,k] = AilM2 
17 



where I = (7,M) is one of the Casimir operators on e(3). Therefore, an 
integrable case occurs at / = 0. The passage to operator representation (30) 
(at / = 0) yields the commuting pair 

Hd = Sx^Dl - {x^ - 8xy + y^)D^Dy + Sy^D^ 

xy — 1 

+ c + 3){xD^ +yDy) + , 

X — y 

-Kd - 4(§ + 1)Hd = ^x^Dl + A{x^ + xy + y^){xD^ + yDy)D^Dy + 4y3z)3 
+ (f + 3)(4(x - yfD^^Dy - (c + 3)(xZ5, + yDy)) 
+ ((x + y)((x2 - l)D^ + (y2 _ i)Dy) - (I + 5)xy + f + 3) . 

4 Spectra 

The quantization in terms of the generators of Lie algebra is universal, but 
the setting of a boundary value problem and computation of the spectra 
depend on the choice of a concrete representation. As an application, we 
consider here the eigenvalue problem for the Euler top on so(3) using rep- 
resentation (6). Recall, that quantization of this model was obtained by 
Kramers-Ittmann [30]. In the case of representation (6), it is natural to de- 
fine the spectrum by condition that eigenfunctions are polynomial. This can 
be compared with Komarov-Kuznetsov paper [22] where the spectrum was 
found for the matrix representation, and the recent Grosset-Veselov paper 
[31] where the spectrum was studied for the representation in elliptic coordi- 
nates and it was shown that coefficients of the characteristic polynomial at 
a given level set j = s are expressed through the so-called elliptic Bernoulli 
polynomials. The spectral problem for the tops on so(4) is rather compli- 
cated and we restrict ourselves by derivation of equations for eigenfunctions 
in the Schottky-Manakov case. 

4.1 Matrix representation of so(3) 

Let us introduce the wave function \m,j) = x^~"^, then representation (6) 
has the following matrix elements in this basis: 

Ui±iU2 = U±, U3\m,j)=m\mJ), 
U+\m,j) = {j -m)\m+l,j), U_\m,j) = {j + m)\m-l,j). (36) 

In this case the condition that j — m is integer follows from the condition 
that the basis functions x^~"^ must be single- valued. If the problem admits 
the time-reversal symmetry m — t- —m then j + m should be integer as well, 
and this implies that j and m are simultaneously integer or half-integer. The 
form of the matrix elements implies that m = —j, ■ ■ ■ ,j, and the condition 
J > follows from the condition that the basis function must be analytic. 
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It is clear from its form that wave function on the orbit j = const is a 
polynomial of degree 2j. 

The wave function and the spectral problem on so(3) can be written as 
follows: 

j 

|^)= C{m,j)\m,j), H\i;) = X\i;) 

m=-j 

Analogously, the wave function and the spectral problem on so(4) read: 
IV') = X] X] C{mi,m2;ji,j2)\mi,m2,ji,j2), H = X . 

mi=~ji mi=-j2 

In the differential representation, the basis wave function on so(4) is chosen 
as |?7ii,m2,ji, J2) = x^'^~"^^y^^~"^^ . Then the general wave function is a 
polynomial of degree 2ji with respect to x and degree 2j2 with respect to y. 

4.2 Spectrum of the Euler top on so(3) 
The eigenvalue problem is of the form 

Hil^f (x) = XtP^ (x) , Cil^^ {x) = j {j + {x) 

where 

H = aiUl + a2Ul + a^Ul C = Ul + + Ul 

We will assume that 01 + 02 + 03 = 0, up to an unessential shift of the 
spectrum. 

Let us use the differential operators representation for so(3) algebra (36) 
and represent wave function in the form 

2j 

i^^ix) = Y,C^{k)xK 
fc=o 

Then the eigenvalue problem is rewritten in the form of recurrent relation 

-^{2j + l-k){2j + 2-k)Cf{k-2) 
+ {\{3{3 + l)-^{j-hf)i-k)cf{k) 
+ ^{k + l){k + 2)Cf{k + 2) = Q 

where f = Si+Sa ^ _ — A — qK _ qX -^j^ boundary conditions on 

^ 01—02' 01—02' J J' 
the left end 

C/(-2) = C/(-l) = 0, C^iO) = = 1. 
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The problem is splitting for odd and even polynomials. 

If j is integer then the boundary conditions on the right end is Cj^(2j + 
2) = for the even polynomials and Cj^(2j + 1) = for the odd ones. 

If j is half-integer then, vice-versa, the boundary conditions on the right 
end is Cj^(2j + 1) = for the even polynomials and Cj^(2j + 2) = for the 
odd ones. 

As a result, the eigenvalues of the Hamiltonian H are zeroes of the 
polynomial 

PjiA) = Cf{2j + l)Cf{2j + 2), degP,(A) = 2j + 1. 

Let us explicitly write down several polynomials -Pj (A) normalized by the 
condition that the coefficient of the leading term A^-'''"^ is unit: 

Po(A) = A, 

A/2(A)=A2, 

Pi (A) = ^(A - 0(2A + e + 1)(2A + ^ - 1), 
P3/2(A) = 3^(4A-9e2-3)2, 

P2(A) = ^(A + 3e)(2A - 3^ + 3)(2A - 3^ - 3)(A - 9^^ - 3)^, 

P5/2(A) = (A3 - 7A(3e2 + 1) + 20^(^2 - l)f. 

All polynomials Pj{^) with half-integer j are full squares in virtue of the 
Kramers theorem [28, p. 225] about the double degeneration of the systems 
with half-integer value of the spin. 

4.3 Matrix representation for the Schottky— Manakov top 

In the Schottky-Manakov case we have two consistent eigenvalue problems 
Hil) = \ip, Kip = The wave function is of the form 

fc=0 1=0 

In the lattice representation, the 5-point equation appears for the eigenvalues 
/x: 

(ai - a2)(2ji + 1 - k){2j2 + 1 - l)C^l%{k -!,/-!) 
+ (ai - a2){k + + l)C^i^j,{k + 1, / + 1) 
+ 2{2{3i-k){j2-l)a^-li)C^l^j^{k,l) 
+ (ai + a2)(2ii + l-k){l + l)C^l^j^{k - 1, / + 1) 
+ (ai + a2){k + l)(2j2 + 1 - l)C^i^j^{k + 1, / - 1) = 0, 
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and 9-point one for the eigenvalues A: 
^al-al){k + l){k + 2)C^l^^^{k + 2^) 
+ {al-al){l + m + 2)C^l^^^{k,l + 2) 
+ 2a3(«2 - ai){k + !)(/ + l)C^l^j,{k + 1, / + 1) 
+ 2a3(Q2 + ai){k + l)(2i2 + 1 - l)C^;'^j,{^ + 1, ^ - 1) 
+ ((ai - a2)2((2ii + l-kf + (2j2 + 1-0' 

- 2(ji + i2 + 1 - k){ji +32 + 1-1) 

- Ul - j2? - jlUl + 1) - J2(j2 + 1)) 

-Aal{{h-kf + {j2-lf) 

+ {a2 + ai)\{ji+j2 -k- If-jiih + 1) - J2(i2 + l)))c^;,^,(fc,0 
+ (al - al)i2ji + 1 - k){2h + 2 - k)C^;^.^{k - 2,1) 
+ {al - al){2h + 1 - /)(2i2 + 2 - l)C^l^,,{k, I - 2) 
+ 2a3(a2 + ai)(2ji + 1 - A:)(/ + l)C^l'^j,{k -1,1 + 1) 
+ 2a3(a2 - ai)(2ji + 1 - A:)(2i2 + 1 - l)C^l^j,{k - 1, / - 1) = 0. 

The boundary conditions are: C^^^-^{k,l) = if the pair (A;,/) lies outside 
the rectangle with the vertices (0,0), (ji,0), (ii, J^), (0, j2)- 

Solutions C^l^-^{k, I) split into solutions on two sublattices: a solution is 
called "even" if it vanishes at odd k + 1, and it is called "odd" if it vanishes 
at even k + I. 

Actually, one can avoid solving 9-point equation: it is sufficient to deter- 
mine the wave functions from the 5-point equation and then the substitution 
into the 9-point one allows to determine the relation between Aj and ^j. A 
plausible answer is that the pairs {Xi,fii) lie on a certain algebraic curve. 

Consider the case ji = j, j2 = ^ as an example. The wave function is a 
superposition of odd and even ones. The even wave function is of the form 

^ = Co + Cixy + C2X^ + C-iX^y + ... 

and coefficients satisfy the boundary conditions C_2 = C_i = C2j+i = 
and recurrent relations 

i(ai + (-l)'^a2)(2j + 1 - A:)Cfc_i + (a3(j - A;)(-l)'= - ^)Cfc 

+ i(ai-(-l)'=a2)(A: + l)Cfe+i = 0. 

Up to the constant factors, one finds 

Co = 0, Ci = //-jQ3, C2 = M^-a3/^-a3j(j-l)-^j'(«i-"2)', ••• 
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The odd wave function is of the form 

^P = Boy + Bix + B2X^y + B^x^ + 

the coefficients satisfy the boundary conditions i?_2 = ^-i = ^2j+i = 
and recurrent relations 

(ai - (-l)'=a2)(2j + 1 - k)Bk-i + 2{a3{k - i)(-l)^ - i2)Bk 

+ iai + i-l)^a2)ik + l)Bk+i = 0. 

Up to the constant factors, one finds 

Bo = 0, Bi=fM + ja3, B2= fi'^ + asfi-a^jij -l)-^j{ai-a2)'^, ■■■ 

The eigenvalues are found from equation 

Pj,i/2{f^) = C{2j + l)B{2j + 1) = 0, degP,- i/2(/i) = 2(2j + 1). 

The polynomial Pj.i/2{^^) is a full square if j is integer, in accordance with 
the Kramers theorem, because then j + ^ is half-integer. 
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